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abstract 

In  a  recent  paper  a  proof  was  given  that  for  a  perfect 
crystal  of  hydrogen  atoms,  described  within  a  certain  model, 
the  free  energy  corresponding  to  localized  one-electron  wave- 
functions  was  less  than  that  corresponding  to  spatially 
extended  one-electron  functions  .  That  proof^  however^  depended 
on  the  assumption  that  the  summand  a^  appearing  in  the  partition 
function  for  the  extended  solutions  monotonically  increases  with 
l  for  l  ^  0  .  The  proof  of  this  monotonicity  is  given  here  . 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt .  Col.,  USAF 

Chief,  Lincoln  Laboratory  Project  Office 
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INTRODUCTION 


1  . 

In  a  recent  paper"*"  (to  which  we  shall  refer  as  I)  ,  a  proof 
was  given  that  for  a  perfect  crystal  of  hydrogen  atoms,  described 
within  a  certain  model,  the  free  energy  corresponding  to  local¬ 
ized  one-electron  functions  was  less  than  that  corresponding  to 
spatially  extended  one-electron  functions  [eq.  (13.25),  i  .e  . 
eq.  (3.25)  of  I] .  These  terms  were  defined  in  the  framework  of 

a  new  variational  approximation  in  statistical  mechanics,  the 

2 

thermal  single-determinant  approximation,  and  under  certain 
specified  limitations  on  the  model  discussed  in  detail  in  I .  This 
result  was  a  central  one  in  that  paper,  "*"  and  therefore  rigor  in 
the  proof  was  strived  for.  However,  the  proof  given  (Appendix  A 
of  I)  depended  on  the  assumption  (made  plausible  there)  that  the 
summand  a^  appearing  in  the  partition  function  for  the  extended 
solutions  [eqs.  (IA.9),  (IA.10)]  monotonically  increases  with  i 

for  i  >,  0  .  The  proof  of  this  monotonicity  is  presented  in  this 
report  . 


1.  T.  A.  Kaplan  and  P.  N.  Argyres,  Phys  .  Rev.  B  1,  2457  (1970)  . 

2.  T.  A.  Kaplan,  Bull.  Am.  Phys.  Soc .  1^3,  386  (1968)  and  Solid 
State  Research  Report  No.  DDC-AD672961,  Lincoln  Laboratory, 
M.I.T.  (1968:2)  p.  53. 
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2 .  PROOF  OF  THE  MONOTONICITY  OF 

For  completeness  we  start  with  the  definitions 


a 


i 


% 


a 


£+ 1 


a 


& 


(2.1) 


(2  .2) 


i  is  an  integer.  (N  is  used  here  in  place  of ~T\  appearing  in  I.) 

N 

We  shall  prove  a^,  ^  af  +  i  f°r  ^  ^  ~2  ~  ■*-  •  Thus  putting 


5  •  y  (2.3) 

and 

& n  g^  =  h  (y )  (2  .4) 

we  need  to  show  that  h(y)  ;>  0  for  0  ^  y  ^  ^  in  steps  of  1/N . 

Equation  (IA.12)  gives 


h(y)  =  N  cosh  ^  +  N  Hn  (1  +  tanh  ^  tanh  xy)  +  ^ — - 

1  +  2y  + 


N 

(2.5) 


But  St,  n  cosh  ^  0,  and,  for  a  ^  0, 


£n  (1+a)  ^  a  - 


(2  .6) 


and 


2 


(2  .7) 


in  (1+a)  a 


Thus 


h(y)  ^  N  (tanh  ^  tanh  xy  - 


tanh2  +  in 


l-2y  1 

T+2y  N(l+2y) 


(2.8) 


But 


tanh  a  ^  a  ,  a  ^  0 


(2  .9) 


so  that  (2.8)  gives 

h(y)  ^  N  tanh  ^  tanh  xy  +  in 

Using  (XA.27),  then 


l-2y  x2  1 

l+2y  "2N  N  (l+2y) 


h(y)  >x  x  tanh  xy  -  in 


1 

N 


or 

h  (y)  ^  7  (y)  -  £  =  hQ  (y) 


where 


7  (y)  =  x  tanh  xy  -  in 


2  3  ! 
p  -  x  .  x  .  1 

t  ~  In  +  — 5  +  n 

3N 


l+2y 

TT2y 


(2  .10) 


(2.11) 


(2 .12) 


(2  .13) 

(2.14) 


The  equation  7  (y)  =  0  arises  in  the  thermal  Hartree-Fock  approxi¬ 
mation .  It  is  easy  to  see  graphically  (see  Fig.  1)  and  one  can 

'V, 

show  analytically  that  there  are  at  most  three  roots  y  =  0,  +  y, 
for  | y |  <  1/2;  also  the  condition  for  the  occurrence  of  three 


3 


roots  is  x  >  2 .  So  if  we  could  drop  £,  then  we  would  know  that 

h(y)  >  0  for  all  y  in  the  range  0  to  y  .  Furthermore,  y  turns  out, 

for  x  as  big  as  100,  to  be  extremely  close  to  1/2,  so  close  that 
1  ~ 

N  (^  -  y)  «  1,  so  we  would  then  have  completed  the  proof.  But  we 

have  £  /  0  .  The  functions  x  tanh  xy  and  £n  j^^y  are  shown  quali- 

2 

tatively  in  Fig.  1,  which  is  drawn  for  x  >  4. 

Our  argument  is  as  follows.  We  will  first  show  that  h(y)  >  0 

for  y  =  ^  i  .e .  for  £  =  ^  -  1  (note  that  for  £  -  N/2,  i .e . 

y  =  1/2,  h(^)  =  -  oo)  ;  that  is,  at  the  first  integral  step  (in  £) 

towards  the  left  from  y  =  1/2,  one  has  passed  the  first  crossing 

of  x  tanh  xy  and  in  +6.  It  is  clear  (from  the  figure)  then 

that  h(y)  will  remain  positive  as  y  decreases  until  y  =  y  ,  (small 

and  positive)  is  reached.  We  will  obtain  an  upper  bound  on  yo  of 

2/N.  Finally  we  will  show  that  gQ  and  g^  are  >  1  and  this  will 

N 

conclude  the  proof  that  a^  ^  a£+l  for  0  x<  &  $  ?  ~  1* 

Putting  y  =  1/2  -  1/W  into  (2  .11) ,  we  have 


h(-r  -  jj)  >x  x  tanh  x  (^  -  ^)  -  in 


7  "  N 


i-* 

~r~ 

N 


-  e 


tanh  -  tanh  . 

- I - n_  _  £n  N(1  1)  _  e 

1  -  tanh  tanh  ^ 


^  x  (tanh  -f  -  ^)  -  £n  N  -  £ 


(2 .15) 
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For  the  last  inequality  we  used  1  )>  tanh  tanh  ^  0 , 

tanh  ^  Clearly  then 

1  1  199  2 

h(-J  -  i)  >  0  for  x  >  200  and  N  <  e  ,  but  x  /N  «  1  (2.16) 

the  range  discussed  below  Eq .  (IA.17)  . 

We  now  find  an  upper  bound  to  the  other  positive  root,  y  , 

of  h  (y)  =  0 .  Put 
o 

1+2y  +  e  =  g  (y)  (2  .17) 

so  y  satisfies 

X  tanh  xyo  =  g(yQ)  (2.18) 

and  is,  by  definition,  the  smallest  positive  root.  Now 
g (o)  =  £  >  0,  tanh  xy  =  0  at  y  =  0,  and  both  tanh  xy  and  g (y) 
monotonically  increase  with  y.  Hence  if  we  replace  tanh  xy  by 
t  (y)  <  tanh  xy  and  replace  g(y)  by  g(y)  >  g  (y)  ,  then 

(yu)  =  g  (yu)  (2  .19) 

where 

yu  >  yc  (2.20) 

(The  latter  may  be  seen  very  simply  by  graphical  means.) 

Using  (IA.27),  we  will  choose 

t(y)  =  xy  -  -i  x3y3  (2.21) 
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Also 


£n  (l-2y)  1_2y 


y  >,  o 


(2  .22) 


Proof :  f  (x)  =  - x  +  in  (1-x)  ;  then  f  (o)  =  0 

£'  (x)  =  - >  0  for  x  >  0 

(i-xr  ' 

.'.  f  (x)  0  for  x  >  0 

So,  using  in  addition  in  (l+2y)  ^  2y,  we  can  choose 

g(y)  =  4y  +  £  (2.23) 

Thus  yu  satisfies  xt  (y)  =  g(y),  which  can  be  written 

2  22143244  — 

(x  -4+2£)y  =  (2x  -4)y  +  x  yJ  -  -j  x  y  +  £  =  h(y)  (2.24) 

Again:  the  i.h.s.  increases  monotonically  and  is  zero  at  y  =  0, 

while  h(o)  =  £  >  0,  so  that  increasing  the  r.h.s.,  h(y),  will 
increase  the  smallest  positive  root.  Clearly,  for  y  <  1, 

h(y)  (2x2-4)y2  +  i  x4y2  +  £ 

so  that  yu  (actually  bigger  than  the  y^  satisfying  (2  .24) )  sat¬ 
isfies 

(2x2-4  +  ^  x4)y2  -  (x2-4+2£)y  +  £  =  0  (2.25) 

We  want  the  root  that  — ►  0  as  £  — ►  0:  So 


7 


(2  .26) 


y 


u 


(x2-4+  2£)  (1  -  v"K) 
2  (2x2-4  +  x  ) 


with 

14  2 

4£  (-i  x  +  2x  -4)  2 

T  =  - ^ - - -  <  1  for  x  200,  P  .*  «  1 

(x  -4+2£)Z  W 


(2  .27) 


_  2 

Also  \  1-r  >  l-r  for  0  <  r  <  1,  and  x  -4+2£  >  0  for  our  parameter 
range ,  so  that 


< 


x  -4-1-  2  £ 

XT 


2  (-^  x  +  2  x  -  4 ) 


2£ 

~~2 


1  - 


2 T 

~7  + 


=  y. 


u 


(2  .28) 


So 


u 


2  -  Ny  = 


,  10  ,4£  2x 

1 - 2  +  X  "  3N 

x  x 


1  "  "X 

X 


+ 


2£ 


>  0  for  the  range  of  interest  (2.29) 


1  .e  . 


ny“  <  2 


(2 .30) 


which  is  the  desired  upper  bound  on  y  . 

2  ° 

For  gQ,  we  use  cosh  x  >  1  +  X-  and  ^  <  £n  (1+x)  <  x  for 
x  <  1,  and  find  that 
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1 


x2/2N 

N  ,0n  cosh  ^  -  £n  (1  +  1  +  x2/2N2 

g  =  e  >.  e 


so  that  g  >  1  for  x  >  200  and  x  /N  «  1 
o 

Again 


(1  + 


~Z 

x 


2N 


(2.31) 


N 


.  ,  x,N  ..  ,  .2  x,N  ~2 

9l  =  <cosh  n>  d  +  tanh  n)  n 


-  1 


f  +  2 


.  ,  .  x.N  ...  6. 

>  (cosh  (1  - 


x2/2N 


>  e 


6/N 

1  +  x2^  ' 


=  e 


S  4-  -  6)  -  ¥  +  H>  '/(1  +  ~2  »  <!  -  H> 


>  1  for  x  ^  200,  |j-  «  1 


(Actually,  to  a  very  good  approximation. 


(>  1  for  x  >  2,  which  is  the 
critical  value  in  the  THFA . ) 
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